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SUMMARY 


This  rc{>oft  is  a  inallicniatical  analysis  of  selected  phases  of  the  Army  Pfiysica! 
Supply  Distribution  System  (APSDS),  The  first  part  is  a  stochastic  analysis  of  e<imoy 
delays?  it  provides  a  statistical  adation  betwflen  distance  travelled  and  tinre  of  travel 
and  incorporates  such  parameters  as  frequency  of  stoppage,  length  of  delay,  and  speed, 
llie  bfcorid  part  aiiaiyswvs  the  problem  of  ci>nvoy  attrition  and  tosses  under  various  hos¬ 
tile  envirunments.  Snch  relations  as  the  expected  number  of  units  lost  and  the  expect¬ 
ed  proportion  of  units  surviving  are  obtained  and  expressed  in  terina^oOnUlaimitnfaer  - 
of  uhitsTn^ie  convoy,  travel  time,  and  various  interplaying  combat-intensity  factors. 
Finally,  a  brief  analysis  of  a  particular  terminal  ojreration  is  given;  expressions  for  car¬ 
go  vessel  turnarmmd  times  and  handlirtg  cost  per  container  arc  obtained. 
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MATHEMATICAL  ANALYSIS  OF  CERTAIN  PHASES  OF  THE 


ARMY  PHYSICAL  SUPPLY  DISTRIBUTION  SYSTEM 

L  INTRODUCTION 

I.  Scope.  The  presanH  study  is  a  malhematieai  analysis  of  selected  phases  of  the 
Army  Physical  Supply  Distribution  System  (APSDS)  in  order  to  attain  a  useful  level  of 
understanding  before  the  necessajy  predictions  can  be  made  about  averaif  system  |ier- 
formatice.  In  addition  to  providing  basic  approaches  to  formulate  and  solve  various 
problems  associated  with  the  dynamics  of  the  physical  supply  distribution  system,  the 
study  demonstrates  the  feasibility  of  approaching  the  problem  on  a  quantitative  basis. 

As  a  result,  it  becomes  possible  to  identify  and  correlate  tbe  relevant  interplaying  input 
parameters. 

As  reported  in  a  previous  study,'  the  APSDS  can  be  conceptualized  as  a  network 
system  consisting  of  nodes  and  links  -  the  nodes  correspond  to  terifiinal  points,  and  the 
links  correspond  to  the  movement  of  supply  between  two  consecutive  nodes.  A  logical 
procedure  would  be  to  study  first  each  link  and  node  and  to  combine  eventually  the  re¬ 
sults  into  a  meaningful  framework  to  arrive  at  optional  solutions. 

The  link  study  is  essentially  an  analytis  of  cargo  retardation  and  attrition.  Sec¬ 
tion  II  is  devoted  to  a  statistical  analysis  of  convoy  delays,  and  Section  III  analyzes  the 
problem  of  convoy  attrition  and  losses. 

Tbe  operation  at  the  node.s  is  more  complex  because  it  involves  a  larger  array  of 
variables  and  covers  a  wider  ^ctruni  of  alternative  systems.  However,  the  cost  and 
time  associated  with  any  terminal  operation  arc  two  accepted  important  characteristics 
which  need  to  be  evaluated.  To  this  end.  Section  IV  illustrates  a  method  of  approach¬ 
ing  the  problem  by  analyzing  the  Sea-Land  Container  Terminal  System;  expre.ssions  for 
cargo  vessel  turnaround  limes  and  handling  cost  per  container  are  obtained. 

When  the  role  of  containers  in  the  APSDS  is  evaluated,  a  fundamental  problem 
area  which  cannot  be  overlooked  relates  to  the  extent  of  penetration  of  containers  into 
forward  areas.  On  the  one  hand,  the  utilization  of  larger  coi'.tainers  presents  economic 
advanl.'};i  -  in  the  handling  cost  at  nodes  and  in  the  transport  cost  between  nodes;  on 

*B.  D.  SivazUan.  “A  Preluninwy  Study  LradinK  to  a  MidiematicaJ  Dcruution  of  «n  Optimum  t'utun-  Attny  Physical 
&ip|dy  Distribution  System,"  USAMEKDC,  Fort  Belvoir,  Virginu  ( Auauat 
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i.he  othfir  hand,  large  containers  could  impede  the  mobility  of  the  distribution  sj'stem 
in  for  ward  areas  and  increase  the^ulnerabilw  of  eo?lV6ys  tfaii^ortri^l:!ic  caigo.  TCerc- 
fore,  one  can  legitimately  wonder  whether  an  optitnum  container  dimension  exists 
which  would  minimize  the  total  cost  of  ojjeration  between  two  successive  nodes.  This 
problem  has  not  been  studied  in  the  pre^nt  report;  however,  the  analytic  results  ob¬ 
tained  to  date  provide  the  necessary'  ingredients  to  approach  thb  unportant  problem. 

IL  RETARDED  CONVOY  MOVEMENT 

2.  lotroduetioa.  Durit^  Uie  normal  course  of  Its  procession,  a  convoy  may  be  im¬ 
mobilized  at  various  Ibne  epochs,  and  each  stoppage  contributes  to  the  delay  of  the 

convoy  in  reaching  its  final  destination,  ^uchjetod^iory  might  Ite  crcatedAyjOUiejcit _ 

tHor^afuraTdr  man-induced  phenomena.  As  an  example,  ah  enemy  attack  might  force 
a  convoy  to  slop  or  to  take  evasive  actions;  the  delay  of  the  convoy  is  in  general  a  func¬ 
tion  of  the  lime  length  of  the  attack.  Blockades  resulting  from  bmrier  encounters  are 
another  cause  for  stoppage;  here,  tlie  delay  will  he  the  time  required  to  clear  and  re¬ 
move  the  barriers  before  normal  procession  can  resume. 

Consider  a  convoy  leavir^  at  bme.  f  ~  o,  origin  A  (x  =  o),  with  tl»e  objective  of 
reaching  destination  B  located  at  distansx;  x  -  X.  Let  v  be  the  normal  convoy  speed. 
We  shall  assume  that:  ( 1)  events  inducing  rodvoy  stoppage  occur  randomly  in  time  ac¬ 
cording  to  a  Poi^H  law  with  parameter  X  (X  is  the  aver*^  number  of  convoy  stop¬ 
pages  per  unit  time);  (2)  the  delay  associated  with  each  stopple  is  a  random  variable 
having  an  exponential  distribution  witli  parameter  p  (l/fi  is  the  average  delay  time  per 
stoppage);  and  (3)  the  occurrence  of  each  stoppage  and  the  induced  delays  are  indepen¬ 
dent  events. 

Wc  shall  derive  an  expression  for  the  probability  that  at  time,  t,  the  convoy 
has  covered  at  least  a  given  distance  and  also  an  expression  for  the  probability  that  a 
given  distance  will  be  covered  on  or  before  a  specified  time.  The  analysis  proceeds  in 
two  steps:  ( 1)  a  time-dependent  model  is  developed;  and  (2)  the  reailts  obtained  are 
then  used  to  develop  a  space-time  dependent  model  involvir^  both  time  and  distance 
as  variables. 

At  the  outset,  without  any  sophisticated  analysis,  it  is  possible  to  obtain  an  ex¬ 
pression  for  the  average  convoy  speed,  v,  for  long  travel  time  (note  that  V  <  v).  If 
we  denote  by'  x  any  distance  interval,  then. 


X  +  expected  delay 


V 


2 


X 


Ji  +  i  .  A  ^  * 

V  V  (I 

1'he  decrenienl  ratio  in  speed  is  X  /  (X  +  n).  Tlicse  results  are  no  longer  valid  if  the 

joumev  tortgULoftho  ^oavcivJa  shorts  ll-wilf  Le-shawn  tlialthe  eurrecf  affswersnare 
functions  of  time. 

3.  Depeadent  Atutys^  We  define  the  following  symbuLs: 

_ — =^fob^iiity^attfae75Snvoyls  floppSTaFtime,  t 

Pj  (1)  =  probability  that  the  convoy  i.s  moving  at  time,  t 

Xdt  =  probability  that  tlie  convoy  will  be  slopped  in  the  interval 

(t,i  +  dt) 

p  dt  ~  probabitily  that  the  convoy  will  start  movii^  in  the  time  interval 
(t,  t  +  dt)  given  diat  it  was  at  a  stop  at  time,  t 

Oj(dt)  =  function  of  dt  such  tfiat  lira  Oj(dl)/dt  =  0  (i  =  1,2) 

dl^O 

We  then  have  the  following  relations: 

P„<t  +  dt)  =  P„<tXl-pdt)  +  P,(t)XdH-0,(dt) 

P,(t  +  dt)  =  P^(t)pdt  +  P,(t)(l-Xdt)  +  Oj(dt)  . 

Carrying  tlie  usual  operation  and  passing  to  the  limit  we  obtain 

-j2j-  =  -pPJt)  +  XP,(t)  ,  (1) 

j  dP,<t) 

**'”*  "TIF"  ^  •  (-> 

It  is  evident  that 

P„{t)  +  P,(t)  =  1  t  >  0  .  (3) 

Initially,  assume  tliat  the  convoy  has  just  started  moving;  then,  l’,{0)  =  0  and  Pj((l)  =  I. 
Using  equations  (1)  and  (3),  we  obtain  the  following  differential  equation  for  P„(l): 
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dP  (t) 

Using  the  initial  condition,  we  ^nd  as  the  solution  of  the  above  equation 

■p  <t)  |rre-«’^>w]-  -  ■  -  - 

Hence,  Pj(t)  -  ^"X+TT  •  (5) 


4.  Expected  Convoy  Speed.  We  define  the  expected  instantaneous  convoy  speed  as 
^[v(t) ] ,  tlie  expectatioit  of  v(t).  Now, 

t{v(t>i  =  vP,(t)  =  V  , 

and  in  the  steady  state 

uv(t)i  =  V . 

5.  Spaee-I^BW  Dependent  Amdyaia.  In  a  spaoe-time  dependent  analysis  (Pig.  1), 
two  state  variables  are  introduced,  namely,  distance,  x,  and  time,  t.  The  determination 
of  two  basic  probability  distributions  is  of  particular  interest:  the  first  is  ?(t,x)dx, 
tlie  probability  ttiat  at  time,  t,  the  convoy  will  have  travelled  a  distance  between  x  and 
X  +  dx;  the  secondis  T(t,  x)  dt,  die  prohability  that  when  die  convoy  is  at  a  distance, 

X,  the  total  elapsed  time  lies  between  t  and  t  +  dt. 

a.  Determination  of  P  (t,  x). 

Define: 


P^(t,  x)  dx  =  probability  that  the  convoy  is  stopped  at  time,  t,  and  lies  be¬ 
tween  x  and  X  +  dx  (0  <  x  <  vt) 

Pj  (t,  x)  dx  =  probability  that  the  convoy  is  moving  at  time,  t,  and  lies  be¬ 
tween  X  and  X  +  dx  (0  <  x  <  vt) 

with  the  symbols,  X,  jj,  and  v,  having  the  same  interpretation  as  before.  We  shall  as¬ 
sume  tliat  the  iiormai  convoy  speed,  v,  is  independent  of  time  and  distance  travelled; 
thus  v  =  dx/dt. 


•-X 
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Fig.  1.  Space-lime  diagram  for  convoy  movement  subject  to  random  delay. 


iipon-enttmt«^ttngthe  vstrioaa  piK^airtiistfe-eveftls,  4li«'foiiewiHg4w^ 

lationa  in  P^(t,  x)  and  (t,  x)  are  obtained: 


P^(t  +  dt,x)dx  =  P^{t,x)dx(l-Xdt)(l-Mdt) 

- - -  -  t_ 

and 

P,  (t  +  dt,  x)  dx  »  P^<l,  x)  dx  ( J  -  Xdt>  Mdt  +  P^<t,  x> djtViit  jadV 


<t,xn-  dx)  dx  (b  ^Jcdt>H-  (]^(dt,  dx),- 


-in 


where  lim  0.(dt,  dx)  « 
dt-0  "*■“« — “  ^ 


i  =  1,2,3 ,4. 


Reiationa  (6)  and  (7)  can  be  written,  re^MJctively,  as: 


P^(t  +  dt,x>  -  P^(t,x) 
dt 


X)  +  XP,{t.  X)  + 


03(dt,dx) 

— at — 


and 


P,(t  +  dt,x)-P,(t,x) 
dt 


=  MP„(t,x>-XPj(t,x)- 


aP,  (t,  x)  dx 
3x  dT 


^  Cl,0t,dx) 
dt 


Let  dt  *  0;  then,  dropping  for  convenience  the  functional  symbolism,  we  obtain  for  the 
two  previous  relations 


aP- 


^  =  A(P  -XP  -v?i!l  . 

dt  «  1  ax 


Ia:t  P(t,x)  =  P^<t,x)  ^  P,{t,x). 


<8) 


(9) 

(10) 


Although  any  of  the  three  quantities,  P,  P^,  Pj  ,  can  be  determined  first  and  the  result 
together  with  equations  (8),  (<l),  or  (10)  can  be  used  to  derive  the  two  remaining  func¬ 
tions,  it  is  easier  to  first  evaluate  the  function  P^(t,x). 
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_EUjQauflating 


at 


+  n 


-it 


-X 


or 


(11) 


Clearly,  each  of  P  and  P^  also  satisfies  equation  (U).  From  the  results  of  the  time- 
dependent  analysis,  Pq(x,  t)  satisfies  the  following  boundary  condition 


/*t 

P>t)dx  =  P„(t)=Y^ 


(12) 


Equation  (11)  is  a  linear  partial  differential  equation  of  the  hyperbolic  type.  'lire  special 
boundary  condition  given  by  equation  (12)  does  not  allow  one  to  obtain  a  solution 
using  the  metliod  of  separation  of  variables  or  by  forming  an  appropriate  Green's  func¬ 
tion.  Using  the  transformation 

w  =  X  and  z  =  vt  -  x, 
equation  (11)  can  be  written  in  the  canonical  form 


,  ji  ap„  X  ap„ 

dw3z  ~"v  ■  V  dz 


(13) 


This  equation  is  amilar  to  one  obtained  by  the  author  in  a  study  of  a  periodic  review 
multi-commodity  inventory  system.*  Assume  a  solution  of  the  form 


A  (Ji 

Po(w,z)  =  c~~'"  ^  .  (14) 

JLj 

1  =  0 

Substituting  equation  (14)  in  (13),  we  obtain  the  following  difference  equation  in  A,: 


*B.  D,  SWulian,  “Inventory  Control  of  «  MulU-Prodoct  System  with  Intersctins  Procurentcnl,"  Working  Piper 
(1966X  Operaliong  Resetrd)  Croup,  Case  Instilute  of  Technology,  Clerdsnd,  Ohio. 
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r^A 


whose  solution  is 


A.  = 


(htY  _J 

\V^  /  C^' 


(r!)2  K 


(15> 


with  being  an  arbitrary  conatant  whose  value  can  be  obtained  by  using  the  boundary 
eonditicn,jBquatioRXl-2).  SubstUutifl^equatioa  ^i5)  itt(14)  we  obtain 


=  A.e-i-'*-"'” 


t  a  0 


-(M 


(2/?“) 


2"  (rt)^ 


Thus, 


1*0  (w,z)  =  A^ev  X^i  wa^ 


where  !„(•)  is  the  modified  Bessel  function  of  orifer  «ero.  The  expression  for  P^(x,t) 


IS 


P„(x,t)  =  I,[|  /^ir(7r-xy| 

To  deteitnine  the  quantity,  A^, ,  we  have  to  evaluate  the  quantity 

P„{t)  =  A^J  1^  [ry  X#i*(vt-x)]  dx. 


(16) 


Using  the  power  series  form  for  !„(•),  we  obtain 

P„(t) 


o  r  =  0 


Since  the  infinite  series  is  uniformly  convergent,  the  interchange  of  the  order  of  ^mma- 
tiun  and  inte^ation  becomes  valid;  hence 
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a  direfiMi^MwofiLwe  obtain  . 


Cojajjaring  this  last  csp^e8l^bR  with  atjuation  C 12).  we  se  ibal  Aj,  -  Ji/v.  faertce 


- — 

We  next  determine  the  quaiUty  P,  (t,  x) .  Freni  equation  (8)  we  have 


Xe»“  R 


Hence 


P,(t,x)  =-ie-^'^(eP‘P,(t,x)) 


1  *  «}  g 

*--6“  »  *  *  ^ 
V  4t 


I.  [2  . 


A  closed-form  expression  can  be  obtained  in  terma  of  (-),  the  modified  Bessei  func¬ 
tion  of  order  1: 


^1(1.5^)  =  ve  » 


V  ®  V  ‘TT^X 


Xx  p(vt»x) 


Findly,  using  equation  (10),  we  obtain  for  P{t,x) 

,  Xx  Mtvt  -  X)  (  f  r-. 


P(U)-;^e-v  V  XIJ2 


Xx  M(vt-x) 


0<X<.1. 

V  ''  J  ) 


In  order  to  determine  P  (t,x),  we  note  that  tlie  probability  that  there  is  no  convoy 
stopple  over  (0,t)  is  simply  e'^‘.  This  is  also  the  probability  that  at  time,  t,  the 
exact  distance  x  =  vt  will  be  travelled.  Thus,  P  (t,x)  has  a  mixed  distribution  with  a 
finite  probability  concentration  at  x  =  vt  and  a  denaty  function  over  0  <  x  <  vt  de¬ 
fined  by  P(t,x).  We  can  thus  write 
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0< x<  vt 


1- 

6 


x~  vt 


x>  vt 


<20) 


1)16  pfobabttity  that  at  time,  t,  a  distance  of  at  least  x  be  covered  Is 
* 

dP{tpj) 


whose  c<Hif%urstion  Is  depicted  in  Fig.  2. 

b.  Detarntitialion  of  7  (Ux)  —  The  Dual  Froi^m.  In  die  study  of  die  dual 
problem,  the  quantity  of  interest  is  F(t,x)  dt  —  tite  probabUity  that  when  the  convoy 
is  at  a  distance,  x,  the  total  elapsed  time  lies  between  t  and  t  +  dt.  The  usi^  of 
the  word  “dual”  will  become  apparent.  It  is  clear  that  the  structure  of  the  function 
a(t,x)  is  such  tiiat 


0 

r*- 

ti 

o 

eH 

1 

V 

(21) 

w(t,x) 

t>* 

V 

since  the  elapsed  time  will  be  t  if  the  convoy  travels  a  distance,  x,  at  the  normal 

V  3C 

speed  without  any  stoppage,  and  the  probabUity  of  mch  an  event  is  e  v  ,  'ffie  func- 
tkm  a(t,x)  is  to  be  determined  from  probainlistic  contideratioos. 

Detine  dt  =  tlie  probability  that  when  the  convoy  is  at  a  distance,  x,  it 

is  stopped  and  the  total  elapsed  time  lbs  between  t  and 
t  +  dt 

and  Sj  (t,x)  dt  =  the  probability  that  when  the  convoy  is  at  a  distance,  x,  it  is  mov¬ 
ing  and  the  total  elapsed  time  lies  between  t  and  t  4-  dt  (t  >  -^). 

It  is  clear  that 

ir(t,x)  =  »^(t,x)  +  fl,(t,x). 

Basic  probabibtic  arguments  yield  the  two  following  relations  in  «^(t,x)  and  ir,  (t,x): 
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?f^(t,x)  dt  =  *)  “  /^dt)dt  +  3Tj(t  •  dt.  X  -  dx)Xdl  •  dt  +  Oj(d(,  dx) 


aitd 


!r,(t,x  +  dx)  dt  =  ffiCt  -  dt,  x}  (t  -  Xdl)  dt  +  rr^(l  -  dt,  x)/jdt  +  Oj((Jt,  dx) 


where 


iiin  o.  (dt,  dx) 


dx-*  0 


“ST 


(}  -  1,2) 


Pagsin^  to  the  limil^ddtapptngdlte  ftinetjoirdlTtotatian,  we  obTain 


3# 


and 


~sr~  “‘^”0  ^*1  ’ 


3^,  3ff 

y -  +  -  _  Xir,  +/iff  . 

3x  3t  *  “ 


nie^j  equations  arc  similar  to  equations  (8)  and  (9);  hence  tr  {t,x)  'satisfies  l)»e  follow¬ 
ing  equation: 


3^ff  .  1  3*e  ,  X  +  o  3ff  .  3ff  _  « 


(22) 


TTie  boundary  condition  for  tins  partial  differential  equation  is  given  by  tlic  normaliz¬ 
ing  condition 


/. 


a(t,x)  dt  ==  1  -  e*^  V 


(23) 


Note  that  P(t,x)  and  ir(t,x)  satisfy  the  same  equation,  hence  the  mtroduetion  of  the 
notion  of  duality.  On  tlie  basis  of  Ute  results  obtained  in  solving  for  P(t,x),  one  may 
assume  that  equation  (22)  has  a  solution  of  the  form: 


^  X*  lJt(yt  -  *) 

—  A  V"  * - 4/'^’“ 


flr(t,x)  =  Ai  e' 


M(vt  -  x) 


..1  h..  -«> 

+  B  •  —  e  V  V 

V 


aIi  [•>  n(vt -'x) 


(24) 
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where  A  and  B  are  arbitrary  constants  to  be  determined  from  tfie  boundary  condi¬ 
tion  of  equation  f23).  Since  the  R.H.S.  term  in  equation  (24)  needs  to  be  integrated 
with  respect  to  t  over  the  range  -^  <  t  <  »  ,  we  evaluate  first  the  foUuwing  quantity: 


/■ 

X/V 


Xx  lt(vt  - 

C”  ~  ■  V -  1 


Let  2  vt  -  X,  then  expressing  as  a  power  series  and  interchanging  the  order  of 
integration  and  summation  we  obtain 


- ™  ^  -  - ^  - 


«  /  XxV  /•  « 

f  fiiiY  dz  ==1 

rl  yv  J  fl  \  n 


\  _  Xx 

r  =  O 


We  next  evah’.ate  the  quantity 

I  -  X 

Vj  (k)  -  I  e"  V  9 
X'  X/ 


X/* 


~/i„[2  /ITT 

3t  1  ^  V 


(I(vt  -  x) 


dt 


Integrating  by  parts,  we  obtain 


v,(x)  =  e'v 


- «) 


I  t)  /~Xx 
l„  /  —  • 

I  V  ^ 


M(vt  -  x) 

V 


(25) 


Using  equation  (25),  we  obtain 
V  (x)  =  I  -  e*  V  , 


Using  equations  (24),  (25),  and  (26)  in  (23)  we  obtain 
v(t,x)Qt  =  .-1+^(1 -e—]= 


(26) 
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tiehce,  A  =  0  and  B  =  v  and  relation  (24)  can  then  he  written  as; 


=  e  V 


.hi.  a 


at 


Xx  #i(vt-  x) 

V  *  V 


I) 


X*  (t{H  ~  *) 
e'  *  ' 


■/ 


[  /Xx  M(vt"-  xT~' 

Xx  ^n<vt  -  x^ 


t  > 


t27) 


'rhis  last  relation  specifies  completety  the  distribution  function  t  (t,x)  as  defined  in 
equation  (21). 

6.  Ihr<d>aiulity  of  Convoy  Reaching  INsiliation  <m  or  bof<we  a  Given  l^ne,  T. 

We  are  now  in  a  position  to  evaluate  Ote  probability,  M(T),  that  the  convoy  will  reach 
its  destination  on  or  before  a  given  time,  T.  Let  x  measure  the.  distance  from  origin 
A  to  destination  B. 


Tlien 


M<T)  =y^djr(t,x) 


More  explicitly,  using  the  results  of  the  previous  anaiysb. 


M(T) 


0 

.Xx, 
e"  V 

Xx 

e”  + 


T<- 


'J*  —  ^ 

V 


ff(t,x)dt 


(28) 


The  configuration  of  the  function  M(T)  is  as  shown  in  Fig.  3.  The  expected  time  to 
travel  a  distance,  X,  is  given  by 


jr  * 

'  Tff(t,x)dt  . 

x/v 


7.  Generalized  Model  for  Convoy  Retardation.  It  is  po.ssible  to  fonnuiatc  a  gen¬ 
eral  model  for  retarded  convoy  movement  in  which  the  parameters  X,'j,and  V  are 
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functions  of  lime,  t,  and  the  distance  travelled,  x.  One  can  eastlyjemonsl  rateLfsii  ^ 
exani|>le*iljalP-^t,x)aiMl  satisfy  equairohs^S)  and  (9)  with  X,#i,and  v 

each  a  function  of  x  and  t.  It  can  dso  be  ^lovnt  that  F^(t,x),  for  example,  satisfies 
the  following  hyperbolic  partial  differential  equation: 


where  f,(t,x)  and  fj{t,x)  are  Itriown  functions  of  X,m,  and  v. 

8.  An  Application  to  a  Prr^lem  of  Equipment  Failure  and  Repair.  The  nmthe- _ 

inatical  mqddsjtliscuascd4hus  far^find  arj  interesting  area  of  ajfpfieation  in  problems  of 
failure  and  repair  of  a  piece  of  equipment.  In  this  ease,  the  interplaying  variables  have 
the  foUowing  interpretation: 

t  =  time 

X  -  number  of  hours  of  operations  logged  by  the  equqnnent 

X  ^  average  number  of  failures  per  unit  time  (failures  are  assumed  to  oc¬ 

cur  m  a  Poisson  fa^ion) 

1/ju  =  average  time  of  repair  (the  repair  time  is  assumed  to  follow  an  expo¬ 
nential  distribution) 


A  specific  application  to  this  problem  relates  to  the  impact  on  cargo  ship  turn¬ 
around  time  of  breakdowns  of  dock  cranes  during  loading  (unloading)  operations.  Fur¬ 
ther  aspects  of  this  problem  will  be  discussed  in  Section  iV, 

Using  Uie  definition  of  M(T),  an  attempt  is  made  in  Appendix  A  to  characterize 
quantitatively  the  concept  of  mobility'  of  an  automatic  equipment  subject  to  break 
downs  and  repairs. 

HI.  CONVOY  ATTRITION  AND  LOSSES 


9.  Introduction.  Cargo  shipments  across  a  link  between  two  given  nodes  take  the 
form  of  convoys  consisting  of  one  or  more  vehicles  carrying  supplies  from  a  given  point 
of  origin  to  a  given  destination.  When  the  convoy  system  in  a  warfare  environment  is 
analyzed,  it  should  normally  be  expected  tliat,  often,  not  all  units  of  the  convoy  (and 
for  that  matter  the  convoy  as  a  whole)  will  reach  their  destination:  in  general,  enemy 
actions  will  tend  to  inhibit  the  normal  supply  operation  and  thus  create  delays,  casual¬ 
ties,  and  losses. 
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In  a  generalized  analyaia  of  convoy  movement,  it  is  necessary  to  study  the  eatu- 

ation  as  a  mised-engagement  problem  involving  four  variables:  time,  t;  the  number _ 

of  units  in  the  convoy,  n;  tlie  number  of  protective  units  escorting  the  convoy,  k;  and 
the  number  of  attacking  units,  b.^  fiy  makii^  suitable  assumptions  about  the  exchange 
rate  between  engaging  forces,  a  Lattcliester>type  equation  can  be  formulated  for  the 
probability  that  at  a  given  time,  t,  the  three  state  vambies,  n,  k,  and  b,  will  take 
on  particular  values. 

'file  eomplexit}'  of  this  problem  is  increased  by  the  numbm-  of  interpiaying 
state  variables  wltich  makes  doubtful  the  derivation  of  amdytie  results^  'ilie  difficutly 
can  be  pardy  overcome  by  comdderii^  die  number  of  cai^  units  to  Ibe  die  only  dgnift- 
cant  vitriabk  and  by  assuming  that  the  convoy  proceeds  in  a  dynamic  environment  (ex- 
lireSible^^oiniTfuhctroirof  time)  causing  attrition.  To  diat  end,  the  environment 
may  be  thought  of  as  causing  diree  types  of  losses; 

a.  Type  1  attrition  induced  by  elemental  attacks  whose  intensity  is  taken  to 
be  proportional  to  the  number  of  units  in  the  convoy.  Thus,  if  there  are  n  units  in 
die  convoy  at  time,  t,  die  probability  that  a  single  unit  will  be  lost  in  the  time  interval 
(t,  t  dt)  is  nX|{t)  dt,  (t)  being  the  intensity  of  attrition  of  Type  1  at  time,  t.  In 
pn'ctice,  a  Type  1  attrition  is  usually  caused  by  the  lack  of  or  the  ineffectiveness  of  con¬ 
voy  protection. 

b.  Type  2  attrition  induced  by  elemental  attacks  whose  intensity  is  uidepen- 
dent  of  the  number  of  units  in  the  convoy.  Here,  the  probability  tliat  a  single  unit  will 
be  lost  in  the  time  interval  (t,  t  +  dt)  is  Xj(l)dt,  Xj(t)  being  the  intensity  of  attrition  of 
Type  2  at  time,  t.  Type  2  attrition  results  if  convoys  are  effectively  protected  or  if  un¬ 
detected  mines  are  tri^ered.* 

c.  Total  annihilation  resulting  from  an  environment  affecting  the  totality  of 

the  convoy.  We  shall  let  A((t)dt  be  the  probability  that  in  the  time  interval  (t,  t  +  dt) 
the  convoy  is  annihilatetl  -  .  . 

Three  typical  situations  will  be  modeled.  Model  1  assumes  that  convoy  losses 
are  due  to  Type  1  attrition  and  an.'ihilation.  Model  2  assumes  that  losses  result  from 
T)'pe  2  attrition  and  annihilation.  Finally,  in  Model  3,  losses  result  from  Type  I  and 
Type  2  attritions  and  total  annihilation.  For  each  of  these  models,  an  expression  for 
P(n,t),  the  probability  that  at  a  given  time,  t,  the  number  of  units  in  the  convoy  is  n, 
is  derived.  Using  this  expression,  it  becomes  possible  to  assess  the  impact  of  the  original 
convoy  size,  N,  upon  such  quantities  as  the  expected  number  of  units  reaching  destina¬ 
tion  and  the  percentage  amount  of  cargo  lost. 

M  Morse  and  G.  C.  KimbaU,  Mfthadi  of  Opemtiom  Research.  J.  Wiley,  New  York.  1951. 

*lbid. 
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a.  Mode!  1  —  Type  1  Attrition  and  Annihilation.  The  equatbm  for  P(n,t)  are 
i*(N,  t  +  dt)  ^  P{N,  t)  ( 1  -  N  X,  (t)dtj  I !  -  ;j(t)dt]  n  =  N 
P(n,t  +  dt)  =  P(n,t)  |1  “nXj(t)dl)  |1 -n(t)dt} 


+  l»(n+  l,t)(n+  l)Xj<l>dt|l 


n  =  J. , ,  N  -  1 


I»<0,  t  +  dl)  =  P(0,0  +  P<l.t)  X,  (l)dt  1 1  -  #i(l)dt]  +  { 1  -  P(0,l)J  fi(t)  dt  n  ==  0 


witfi  the  initial  conditions  P(N,  0)  =  I  and  P(n,  0)  =  0 


n  =  0, . .  . ,  N  -  1 


The  above  equations  reduce  to  the  following  system  of  differential-difference 
equations; 


dP(N,t) 


=  -  [;i(t)  +  NX,(t)}  P<N,t) 


=  -  I;i(t)  +  n  X,  (t)l  P<n,t)  +  (n  +  1)  X,(t)  P(n  +  l,t)  N  =  I, . . .  N  - 1 

^  iu(t)  [  1  -  P (0,t)]  +  X,  <t)  P(l,t)  n  =  0 

These  equations  can  be  solved  recursively  to  yield 

f  1,2 . N 

W,.)  =  I-e'-t"'"”'  [l-.'-t'*-'*''’"]''  „  =  0 

The  expected  number  of  units  reaching  destination  is 


which  is  proportional  to  the  number  of  units  originally  in  the  convoy.  The  expected 
number  of  convoy  units  lost  is 


./‘iX,  (d)  +  wfl)Ndl 

N-  f(N)  =  n[i  -  e*t  ‘  1. 
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^id  the  expected  proportion  loss  in  cargo  ia 


|{N)  I  -  e  ’ 

which  is  independent  of  the  number  of  units  in  the  convoy. 

. _b.  -  Model  2^  Type-2;4ttritwn  and  AnnStUili^-  The-e^atioiis-for-l<(«it)  aac: 

P(N,t  4  dt)  =  P(N.t)  n  -  Xj<t)  dtn I  -  ;i{t)  dt  I  n  =  N 

P<n,  I  +  dt)  =  P(n,t)  I  i  -  X, <t>dt]  ( 1  -  m  dt) 


+  P<n  4  1,  t)  Xa  (t)dt  {1  -  M(t}dtl 

P(0,t  +  dt)  =  P<0,t)  +  P(l,t)Xj(t)dt)l -M(t)dt) 
+  [l-P(0,t)|  /i(t)dt 


n  =  1, . . .  ,N  -  I 


n=  0 


with  the  initial  conditions  P(N,0)  =  I  and  P(n,0)  =  0  m  =  0,  1  . , . ,  N  - 1.  This  reduces 
to  tlie  following  system  of  differential  difference  equations: 


=-[X,(t)4,.(t)I  P(N,t) 


+  P(n,t)4X,(t)P(n4l,t) 

=-u{l)P(0,t)  +  M(t)  +  Xa(t)P(l.l) 

These  etjuations  may  bo  solved  recursively  to  yield: 


Uit  iNwi  ^ 

^  X,  {.0)6.9  iXa  (9>  +  ih9)  1  id 


(N  -  n)! 


M  =  N 

P(0,t)  =  1  -  y]  P(n,t) 

JSmmm 

M  =  l 

0(«)  =  £ 


n  =  N 

n  =  1,2, . . . ,  N  -  1 
n  =  0 


n  =  l,2 . N 


n  =  0 
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Hie  expected  number  of  umts  reaching  de^jination  jfL 
0=8 


ft  «  o 


- - - - ‘ 

J  "  {N  -  n)! 

ft  *“  i 

•  G{N  - 

O 

In  Appendix  0,  we  sliow  that  the  K.H.S.  term  is  an  increasing  and  convex  function  of  N. 
The  expected  proportion  loss  in  cargo  is 

g(N, .  o(N-i)]  . 

We  ^ow  in  Appendix  B  that  this  is  a  decreasing  function  of  N. 

e.  MocM  3  —  Types  1  and  2  Attritson  and  Aniiihilati<m.  The  equations  for 
P(n,t)  are 

P{N,t  +  dt)  =  P(N,t)  1 1  ♦  N  X,  (t)dt]  1 1  -  Xj  (t)dtj  [  1  -  ^«{t)dt]  n  =  N 

P<n,t  +  dt)=  P(n.l)  Il-nX^{t)dt]  {1-X,{t)dt]  [l-t<(t)dtj 

+  P(n+  l,t){n+  i)X,(t)dtl  {I  -X,(t)dt}  [1  -p(t)dtl 
+  P<n  +  l,t)  [  1  -  (n  +  1)X,  (t)dt}  X,(t)dt{  1  -  #r(t)dt]  »  =  1,2,...,  N  -  1 
P{0,t  +  dt)  =  P{0,t)  +  P{l,t)  X,  {t)dl  ( 1  -  a(t)dtl  [  1  -  X, (t)dt) 

+  P(l,t)  f  1  -  X,  (t)dtl  XjWdtfl  -  i«(t)dtj  +  [I  -P{0,t)J/s(l)dt  n  =  0 

"nie  initial  conditions  are  P{r'l,0)  =  1  and  P(n,0)  =0  n  =  1 ,2,...,  IN  -  1. 

The  corresponding  differential-difference  equations  are: 

+  X,(t)  +  NX,  (t)]  P(N,t)  n  =  N 
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^  .  {1)1  - 

+  (Xj(t)  +  (n+l)Xj(t)]  P(n  +  l.t) 

=  -ii{t)  P(OJ)-Mt)+  jX,(t)  +  54(t)]  P{1,0 


n  =  I,  2, . , . ,  N  -  I 

ft  ~  0 


Let  Vi 


ri 


<?j<t)  =  J*  da 

O 

Then,  we  obtain  recuraivety 

P(N,t)  =  ^  NX, {9)5  <10 

ft/w  1  .V  -i*l«0)+^(0)*(N-»)X  (0)Jdfl 

P{N  -  l,t)  =  e  **»  ^  ‘  (v»j(t)+ N^,(t)] 

P(N-2,t)  =  „-/.l«»). 


/ls<“)  + 


(N-l)X,(u)l  [.^j(u)  +  NvJ,(u)}du 


Unfortunately,  closed-form  expressions  cannot  be  obtained  for  P(n,t).  For  the  case 
when  X,  (t),  X}<t),  and  p(t)  are  independent  of  t,  a  simpler  expre^on  can  be  denved. 
Let  a  =X^/X,,  then 


P(0,t)  =  1  .  + 


r(a  +  i).r(N+  1)  ®  J  ^  “iV  ®  /  J 


The  expected  number  of  units  in  the  convoy  is 

n  ==  N 

-(Xj  *  #1)>  P(^ 


E 


"X 


r(n  +  a 


a+lKN-n)!  \®  /  ^  ®  } 
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Apraafiaxieiiionsti^that- N  andtiiat  g<N),  the 
expected  propertion  in  cargo,  k  a  decreasing  Unction  of  N  eiudeg  ua  for  the 
present. 

10.  Coat  CoiakiefatioB  in  tiie  Moventent  of  a  Convo)'.  1'he  three  basic  cost  elements 
that  must  he  accounted  for  in  the  of  cai^o  are: 

a.  Ihe  cost  of  transpc^tion  winch  is  an  increasing  function  of  the  number 
of  units  in  the  convoy,  N. 

_ I'he  coat  of  providing  a  protective  CacorLio  thttconvay:  in^generai^il  kes- — 

pected  that  this  cost  be  made  up  of  two  components,  a  fixed  component  and  a  variable 
component  which  are  increasing  functions  of  N. 

c.  The  cost  assoeiated  with  the  loss  of  cargo  units;  this,  again,  is  expected  to 
be,  in  general,  an  tncreaang  function  of  N. 

The  conf^ration  of  the  first  two  cost  functions  has  to  be  determined  em¬ 
pirically;  aldiough,  as  a  fir^  good  approximation  they  can  be  assumed  to  be  linear  in  N. 
The  third  cr»t  function  is  easily  determined  hom  previous  analyses  since  it  will  be  equal 
to  a  constant  dme-the  rpiantity  [N  •  f(N)} .  It  then  becennes  possible  to  expre^  the 
totsd  cost  of  moving  a  convoy  as  a  function  of  the  number  of  units  in  the  convoy.  Finally, 
it  might  be  a  matter  of  interest  to  determine  the  cost  of  moving  a  single  ton  of  commod¬ 
ity  between  given  nodes. 


IV.  ANALYSIS  OF  A  TERMINAL  OPERATION  -  THE  SEA-LAND 
CONTAINER  TERMINAL  SYSTEM 

11.  Introduction,  A  terminal,  or  node,  can  be  defined  as  a  transfer  point  where 
cai^o  materiel  experiences  one  or  more  of  the  following  operations:  handling,  storage, 
unidzalion,  de-unittzadon.  In  general,  the  mathematicai  analysis  at  nodes  will  depend 
upon  the  porhcular  terminal  confipiration  and  operation.  In  this  chapter,  the  Sea-Land 
Container  terminal  operation  is  analyzed. 

Containers  (trailers)  mounted  on  truck  chassis  wait  in  a  marshaliiirg  area.  When 
a  container  vc^l  docks,  a  stack  of  trailers  in  one  of  the  vessel  cells  is  unloaded  one  at  a 
time  onto  waiting  truck  chassis  by  one  or  two  gantry  cranes.  Ttie  chassis  arc  pulled  over 
to  the  marshallit^  area  by  tractor.  Following  this  initial  operation,  tlte  crane  unloads 
an  export  container  in  the  empty  cell;  and,  on  its  return,  tire  crane  lifts  an  import  con¬ 
tainer  from  an  adjacent  ceil  and  deposits  it  on  tlie  truck  chassis  thus  initiating  a  one-on, 
one-off  cycle.  Tractors,  in  the  meantime,  move  back  and  forth  between  the  dock  and 
the  mar^alling  areas  pulling  alternately  import  and  export  container.^. 
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Because  of  the  difference  between  the  operaBng  cycle  Bmes  of  craac  aad_trac. 
ttw,  It  beeomeff  nee&asafy  lb  defermfite  tffe  optimura  itumber  of  toetors  ui  operation. 
Simulation  studies  of  Uiis  problem  were  performed  by  the  National  Farts  Council.^  In 
this  chapter,  we  develop  two  mathematical  models  to  describe  tire  probiem.  In  tire  first 
model,  the  crane  cycle  time  and  the  tractor  cycle  time  arc  a^unted  to  be  determini^ic 
variables;  in  die  second  luodci,  die  same  variobies  are  amimed  to  be  stochastic. 

Four  basic  c^rations  to  consider  are; 

•  the  tractor  travel  time 

•  die  tractor  waiting  time _ 

•  the  lift-on,  lift-off  crane  operation  time 

•  the  crane  waiting  time 

By  minimising  the  total  cost  of  handling  a  unit  container,  the  optimum  number 
of  tractors  to  be  sequenced  with  the  crane  can  be  determined. 

0.  The  Deterministic  Model.  Let: 

n  =  number  of  tractors  aligned  to  each  crane 
gc  -  operation  cycle  time  of  crane 

s’"  =  tractor  travel  time  (assumed  the  same  for  all  tractors) 
w®  =  waiting  time  of  cranes 
~  waiting  time  of  a  tractor 
T  =  total  cycle  time 

=  cost  of  operating  die  crane  per  unit  time 
Cj.  -  cost  of  operating  a  tractor  per  unit  time 

Then  T  =  ni?'  +  w' 

= 

Hence  (n  -  I )  fl®  +  w*  =  . 

It  is  evident  that  if  the  time  parameters  involved  do  not  change  over  time,  then  either 
the  crane  wails  or  the  tractors  wait  These  two  cases  are  illustrated  in  Fig.  4. 


^NtHoiul  PtrU  Council  Ketcarefa  and  Technical  Butktin  Na  21967, 17  North  Audky  Si.,  London,  W,  L 
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TRACTOR  I 


II'I  I  r'lTl'ii’"  •^Tr’*1'~Tili'Vnii‘'ini1i*l‘li‘il‘i*iririMi  Ii*p'li  iil*TiliTlrtfflhiaiin'T1t'iiilViTi  iif'IViTiin’i  nr' ii  n  •  '  ill  I  |-  'ii'i  <  <  '  II  •|*'II| 


(1)  IVoporti<Mt  of  Waiting  Time.  If  w*=0,  then  w*‘#=0,  and  T  =  n^ 

=  eycle  time.  Each  tractor  will  then  wait  an  <■  niount 

=  (n  - 1)5“  -  5^  (n  -  1)5“  >  S’" 

and  the  proportion  of  tractor  idle  time  per  cycle  is 

_  {n-  i)6^  -5^  . 

T  "  n5' 

If  now  w“  0,  then  w^  -  0  and 

T  =  5’^  +5=  .  _  _ _  _ _ _ 

Itie  crane  will  then  wait  an  amount  equal  to 

w*  =  5'*'  -  (n  -  1)5^  5'*'  >  (n  -  1)5“ 

per  cy  cle,  artd  the  proportion  of  crane  waiting  time  is 
w^  _  5“^  -  (n  -  l)lf 

In  Fig.  5,  the  waiting  time  per  cycle  is  plotted  as  a  function  of  the  number  of  tractors 
in  operation. 

(2)  Loading  —  Unloading  rate.  The  loading  —  unloading  rate  is  measured 
by  the  number  of  containers  handled  per  unit  time;  the  expression  for  this  rate  will  de¬ 
pend  on  whether  the  crane  waits  or  the  tractors  wait.  It  is  given  by  the  following; 

gT 

n  <  +  1  (crane  waits) 


n>  I  (tractor  waits) 

6 

Thus,  if  the  crane  is  allowed  to  wait,  the  loading  rate  is  proportional  to  the  num¬ 
ber  of  tractors  involved;  otherwise,  if  the  tractors  have  to  wait,  the  loading  rate  is  con¬ 
stant.  This  is  illustrated  in  Fig.  6. 

(3)  Total  Turnaround  Time  of  Cargo  VesseL  Denote  by  N  the  total  num- 
l»er  of  containers  to  be  loaded  (unloaded);  then,  neglecting  the  initial  first  cel!  clearing 
operation,  the  total  turnaround  time  is  given  by 
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1 


/  n(e^  +Y) 

I  " 

KM 


if  n  <  1  + 


6^ 


if  n  >  1  +- 


aitd  thi$  is  ilisplayed  as  a  function  of  the  number  of  tractors  in  Fig-  7, 

<4>  Nue^ier  of  Tractors  in  Operation.  Conmdering  as  our  ob> 

jeclive  function  Uie  cost  of  handling  a  container,  say  C(n},  then 


Q  +  c^(0T 


C(n) 


if  (n  -  1)  < 


if(n-  1)  > 


Tire  function  C{n)  is  sliown  in  Fig.  8,  and  its  analytic  properties  are  investigated  in  Ap¬ 
pendix  C.  There,  it  is  shown  that  tlie  optimum  value  of  n,  say  n*,  which  minimizes  tire 
function  C(n),  is  such  that  one  of  the  following  conditions  are  satisfied: 

0'*'  C. 

(e=Y  >  ‘“teger). 

(a)  Crane  will  wait: 


^  -(C-0)  +  y(C  -^0)^  +  4C 
<  2 
(b)  Tractors  wait: 


<  n*  <  1+0 


-{C-e)+^/^+8f  +  4C 
(l+0)<  n*  <  1  +— !^ - - 

(c)  No  waiting  for  tractors  or  crane; 
n*  =  1  +  0 

(d)  Double  solution  (either  crane  waits  or  tractor  waits): 

-(C-0)  +  /(C  +  0)^  +  4C” 

“*  ~  2 
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Fig.  8.  Handling  cost  per  container  as  a  function  of  the  number  of  tractors  used 
(deterministic  case). 


It  is  evident  that  full  utilHestlott  of  the  handling  system  in  a  rigidly  si^ueneed 


^ _ _ _ _ _  -  — 

operation  Vs  achfcvedlBy^lieFn^  the  quantity  &  as  an  integer;  then,  under  opti¬ 
mum  operating  conditions,  the  cost  of  handling  per  container. 


i5a:litiear  fia«;tion-of  the  crane  cycie-time'afid'tiieTr^tdr'thiva  time  —  the  optimum 
value  of  n  being  (ff  +  1),  a  quantity  independent  of  the  cost  parameters.  Hm^ever, 
in  generd,  ^ch  scheduilng  catutot  be  adtieved  even  within  a  single  ioading-mtloadisg 
ofM^ration  of  a  vessel  faeomse  tite  quantify  is  a  function  of  the  distance  travelled  by 
^  tiractor  and,  hertce,  a  function  of  tiie  locatimi  coordinate  of  tim  iftar^_ 

shailingama. - 

b.  The  Stochastic  Model  In  practice,  a  rigidly  sequenced  operation  is  impos¬ 
sible  to  achieve.  The  tractor  travel  time  will  vary  between  tractors  and  between  travels 
vdille  die  cycle  time  (qmration  of  die  crane  will  have  inherent  delays  due  to  hatch  re¬ 
moval,  hat^  replacement,  latching  operations,  and  movement  of  the  cranes.  We  tiiall 
asaune  that  arrival  of  the  tractors  to  discharge  berth  is  Poisson  distributed  with  inten¬ 
sity  of  arrival  X;  while  die  time  for  the  crane  to  complete  the  one-on,  one-off  opera¬ 
tion  has  an  exponential  distribution  with  parameter  m  (note  that  0^  =  1/fi  and  0^  ~ 
l/X).  We  have  contidcr^,  here,  an  extreme  case  of  randomness;  and,  althou^  the  as¬ 
sumptions  about  the  statistical  distribution  need  verification,  the  results  obtained  can 
be  used  as  bounds  for  tire  values  of  Urn  variables  involved. 

The  problem  under  study  is  similar  to  the  classicat  repairman  probfem.* 

Let  P(ni,  t)  be  the  probability  diat  at  time,  t,  there  are  m  tractors  waiting  in  line  by 
the  crane  (tractors  being  serviced  on  a  first  come  first  served  basis)  m  <  n,  then 


P(0,t  +  dt)  =  p(0,t)(l-nXdt)  +  P(l,t)i«dt[l-{n- l)Xdt] 

P(ai,t  +  dt)  =  P(m,t)(l  -  ftdt)  [  1  -  (n  -  m)  XdtJ 

+  P(m-l,t)(l-jtidt)[n-(m-  l)Xdt]  +P 
+  P(m  +  l,t)pdt  [1  -  -  (m  +  l^Xdt]  m  =  1,2, . . . ,  n 

These  relations  reduce  to  the  following  system  of  differential-difference  equations: 
^?^  =  -nXP(0,t)  ,rP(l,t) 


m  =  0 


m  =  0 


t'cUer,  Introduction  to  Probability  Theory  and  iu  Applkalioiu,  VoL  1,  J.  VSey,  New  Yoik,  1949. 
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i  i 


ti 


tl 


^  ^mS)  7 

=  (It  -  m  +  1)X P{in  -  l,t)  -  +  (n  .  m)xj  p^n.t)  +  P(m  +  l,t) 

m  =  1,2, . . , ,  n 

lo  Aes^dy state,  Jet  P„  -  P(m,!}.  (m -o,  1.. ..  ,a),  then 
0  -  -nXP^  +  mP, 

0  =  (a>m  +  l)XP,^  j -(^  +  (n.ni)X)P,^ in=l,2,...,i 


m-O 


It  m  easy  to  verily  that 

p  =  fiT  .  »7  p 

®  Kl*  /  (n  -  m)}  o  ’ 


m=  1,2,..  .,11 


and  P  s  - 

o  /..%a  m=a/  ..vm 


=  !>.<») 


feOTdled  per  unit  fame  by  the  crane  while  in  operation  (not  Idie)  is  it  follows  that 
the  expected  number  of  co.itainers  handled  per  umt  time  is  (1  -  P  ) 

(1)  Expected  Tomaround  Tiiae  {Fig.  9).  This  is  given  by 
N 


II -wu 

where  N  is  the  number  of  containers, 

u-  ^  of  tractore  in  operatic  {Fig.  10).  Considering  as 

the  objective  ftincbon  the  expected  toUl  cost  of  operation  per  container,  C(n),  we  have 


C(n)  = 


_  C^+nCj 


^e^optimum  number  of  tractors  can  then  be  determined  by  minimiaing  the  function 
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NUMBER 

OF 

tractors 


Fig.  9.  Tolal  ship  iurnaround  lime  as  a  function  of  the  number  of  tractors 
{stochastic  case). 


5  6 


8 


®  NUMBER 


OF 

TRACTORS 


Fig.  10.  Expected  handling  cost  per  container  as  a  function  of  the  number 
of  tractors  used  (stochastic  case). 


12.  Effect  of  Crttw  Breakdown  on  Ttimaroutid  Hate.  To  this  point,  we  have  as* 
samed  that^  e^ipuetits  thv^ved  in  ffie  loaiiui^-uiiload&tg  dpefation  ^  not  subjeet 
to  breakdowm  the  incorpomtion  of  tlug  factor  and  its  impact  upon  turnaround  tune 
is  a  at^fkant  probtem  that  we  dtidl  now  discus.  Since,  in  a  Sea-Land  terminal  system, 
crane  is  Uie  bottleneck  eijuipmeni,  we  shall  study  the  effect  of  crane  breakdown  on 
turoaroond  time.  Aasuine,  fust,  tiiat  a  m^e  crane  is  operating.  Let  the  capadty  of  the 
vea^be  St  containers,  andlet  -T-  be  die^maaJiiiunt-pemiiseftk^binua^u^tune. 
Under  normal  operation  (no  broiJtdowns),  the  drip  tiimarotmd  tune  wiU  be  dictated  by 
die  number  of  containers  loaded  on  the  drip;  however,  in  case  of  breakdown,  ihe  d^p 
bmunound  time  m^t  well  be  ihctated  by  T  (F^.  1 1).  Let  the  ioading  rate  be  v  eon- 
tainers  per  unit  time,  and  as^une  that  breakdowns  occur  randomly  in  tune  according 
tea  Poisaotrkw  wiiitinleiisity^T^ttepairtBiiebeexj^SSBnfi^y^isTribul^  with 
parameter  iu  .  The  ^tuation  is  similar  to  die  convoy  retardation  problem  discussed  in 
Section  11.  it  is,  thus,  possible  to  obtain  an  expression  for  the  expected  turnaround  time 
as  a  function  of  N,  T,  v,  X,  and  n. 


V.  CONCLUSIONS 

13.  Cotudusknis.  The  mathematical  analysis  so  far  performed  has  been  restricted  to 
selected  phases  of  the  ATSDS  such  as  the  movement  of  convoys  across  links  and  the 
'iteration  at  a  node  having  a  special  configuration.  The  analysis  is  not  claimed  to  be  ex¬ 
haustive;  in  fact,  the  solution  to  some  of  Ae  mathematical  problems  encountered  etudes 
us  at  present.  Nevertheless,  the  results  obtained  so  far  constitute  part  of  the  basic  inputs 
necessary  to  characterize  the  optimal  configuration  of  the  APSDS. 
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APPENDIX  A 


CONCEPT  OF  MCaiUTY  OF  AN  AUTOMOTIVE  EQUIPMENT 
SUBJECT  TO  BREAKDOWN  AND  REPAIRS 

The  mobility  of  an  automotive  piece  of  equipment  operating  under  prescribed  con¬ 
ditions  can  fae  defined  as  the  probability  that  the  equipment  starting  from  a  given  origin, 
G,  and  moving  along  a  giv  nj  path  will  reach  a  predetermined  destination.  A,  on  this 
path  nn  or  befrire  a  specified  time,  T : 


Tile  quantity  M(T)  defined  in  Section  II  by  expression  (28)  can  thus  be  used  to  meas¬ 
ure  mobility.  Note  that  M(T)  incorporates  the  following  four  ba.sic  factors: 

1.  normal  speed  of  motion,  v 

2.  distance  travelled,  x 

3.  frequency  of  breakdown,  X 

4.  average  repair  period,  1 A 

Consider  two  motor  equipments  E,  and  moving  from  O  to  A  along  OA  at 
the  same  speed,  v.  Assume  that  Ej  and  Ej  are  subject  to  the  same  law  of  failures  but 
that  E,  is  supported  by  a  better  repair  system  than  E^;  then,  is  more  mobile  than 
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APFKNDIX  B 


vv(“  have,  for  iar«e  N, 


f(N)  A  +  Ne'^* 
where  A  <  0  is  a  constant. 

Next,  we  study  the  function  g(N); 
g(N)  -  [g{N>.  Mg<N-  1)] 


-  <Xt>^^'  (1N  + 1)G(N- l)-NGfN> 

["irri)!'  -  '-'n{n  +  i) — 


Finaily, 


Xte^^"»‘  f  mf 
N{N+'l)  STT-I^ 


+  G(N  -  1)  -  NG(N  -  1) 


+  NG(N  -  1) 

(Xt)^  \ 

■  0^  -  or  j  ' 


g(N  +  l}-g(N}  =  - 


NfN+iyGCi^  -  1)  <  0  . 


Therefore,  g(N)  is  a  decreasing  function  of  N.  The  convexity  and/or  concavity  of 
g(N )  can  be  established  as  follows: 


!g<N^2)-g(N  +  l))  -  [g(N+l)-g(N)]  =Xte'^^"^^‘ 


G(N)  C{N-1)] 
(N+l)(N+2)  N(N  +  1) 


=  -XI 


-(X  +  t 

““fTTF' 


2G(N  -  1) 


(Xt)'" 
(N-  I)! 


Since  tliis  last  expression  ean  be  either  positive  or  negative  depending  on  the  specific 
values  taken  by  Xt  and  N,  it  is  expected  that,  in  general,  g(N)  will  not  exhibit  any 
particular  convexity  or  concavity  property. 


40 


APPENDIX  C 


ANALYTIC  PROPERTIES  OF  THE  FUNCTION  C(n) 

Herein,  we  develop  the  conditions  which  will  dictate  the  uplimtiin  number  of  tractors 


to  be  used  when  the  quantity  6 


is  not  an  integer  (Case  2). 


It  is  evident  that  the  optimum  number  of  tractors  to  use  is  either  n*  =  [5  + 1 )  or 
n*  =  [0+2]  or  both,  where  jxj  is  a  ay  mboJ  denoting  the  ssaailest  integer  in  the  real 
quantity,  x.  The  appropriate  value  of  n*  can  be  determined  through  considerations  of 
tlie  objective  function  which  is  taken  to  be  the  total  cost  for  handling  per  unit  con¬ 
tainer. 

The  convexity  of  the  objective  function  guarantees  the  existences  of  an  optimum  n*; 
however,  in  this  case,  a  double  optimal  solution  is  possible.  Further  analysis  to  deter¬ 
mine  preferability  of  crane  waiting  rather  than  tractors  (or  vice  versa)  proceeds  as  fol¬ 
lows: 


•  '  1 

n*  rfVi 

Assume  that  at  the  optimal  point,  when  n  =  n*,  the  crane  has  to  wait;  then,  n*  <  0 
+  1.  Clearly,  if  (n*  +1)  tractors  are  used,  the  tractors  will  be  w'aitiiig  and  the  follow¬ 
ing  inequalities  should  then  hold: 

<  (0^+0^)^.^  +  C.,)<  0^[q+(n»+  l)C.j].  <1 

The  L.II.S.  inequality  yields 
n*  <  e  +  1. 


The  R.II.S.  inequality  yields 


(■+»)<  . 
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or 


n*  +  (c  -0)11*^-  c(J  +0)  >  Q  . 


Since  the  corresponding  quadratic  equation  has  two  real  roots  of  opposite  signs,  this 
last  inequaiity  will  hold  for 


n*  > 


-(c-0>  +  y^<c  +  of  +  4e 


2 


It  is  easy  to  verify  tliat  the  R.ElJs.  term  of  this  inequaiity  is  greater  than  0.  Hence,  the 
set  of  inequalities  (I)  can  equivalently  he  written  as: 


+  /  (c  +  of  +  4c 

e  < - X - <  n*  <  1  +  0. 


(2) 


Assume,  now,  that  at  the  optimal  point,  when  n  =  «*,  the  tractors  have  to  wail;  then, 
n*  >  1  +  0  . 


Since,  if  (n*  -  1)  tractors  are  used,  the  crane  will  have  a  waiting  time,  it  follows  that: 

(C,  +  )  <  f  (C,  +  n*C^)  <  (C,  +  (n*  -  l)C^)  .  (3) 

The  L.H.S.  inequality  yields 
n*  >  1  +  0  . 

The  R.II.S.  inequality  yields 

(n*  -  l)(c  +  11*)  <  (1  +  0){c  +  n*  -  1)  , 

+  (c  -  2  -  0 )  n*  -  2c  +  I  -  0c  +  0  <  U  . 


'I'lie  corresponding  quadratic  equation  can  be  shown  to  have  two  real  roots  given  by 


-(c-&)  ±  '+  4); 


(^onw^quently, 

~(c-e)  -  /(e  +  of  +  4c  +  /(c-i  fff  +4~ 

I  +  - - - <  n*  <  I  + - — - - - 

ll  ran  he  verified  lhal  the  left-most  quantity  of  the  above  expression  is  less  than  i  +  ; 
hence,  the  set  of  inequalities  (3)  can  equivalently  be  written  as 

-(c-d)  +  j' (c  +  of  +  4c 

(1  +9)  <  n*  <  I  + 


2 


